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Abstract 



We investigate, using the noise reduction technique, the asymptotic uni- 
versality class of the well-studied nonequilibrium limited mobility atomistic 
solid-on-solid surface growth models introduced by Wolf and Villain (WV) and 
Das Sarma and Tamborenea (DT) in the context of kinetic surface roughening 
in ideal molecular beam epitaxy. We find essentially all the earlier conclu- 
sions regarding the universality class of DT and WV models to be severely 
hampered by slow crossover and extremely long lived transient effects. We 
identify the correct asymptotic universality class(es) which differs from earlier 
conclusions in several instances. 
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Typeset using REVT^ 



Kinetic surface roughening of nonequilibrium growth models, particularly under solid-on- 
solid (SOS) growth conditions, remains a subject of considerable interest and activity in spite 
of a great deal of theoretical and experimental research in the topic during the last decade 
While the generic nonequilibrium surface growth universality class (in the situation 
allowing overhangs and bulk defects, i.e. under generic non-SOS conditions) is theoretically 
accepted to be the KPZ universality (although the asymptotic KPZ universality in specific 
models may be masked by extremely slow crossover effects), there is no such consensus for 
SOS growth models where seemingly innocuous small changes in local growth rules appear 
to lead to different dynamic universality classes. In fact, some of these nonequilibrium 
SOS growth models, introduced in the context of mimicking ideal molecular beam epitaxy 
(MBE), do not even seem to obey self-affine dynamic scaling behavior, instead exhibiting 
nonuniversal anomalous and multi-affine scaling. It may therefore be questioned whether the 
universality class concept is useful in SOS growth models or even applies to nonequilibrium 
SOS growth models. 

In this paper we address this confusing situation plaguing our understanding of the uni- 
versality class of discrete SOS nonequilibrium growth models by concentrating on a specific 
class of surface diffusion driven stochastic growth models which mimic in a drastically sim- 
plified manner, low temperature molecular beam epitaxial growth. This class of models has 
been referred to as 'limited mobility' nonequilibrium growth models because the deposited 
atoms are only allowed to diffuse (obeying certain specific local diffusion rules) at incidence, 
and all the other atoms in the growing film, except for the most recently incident atom, 
do not diffuse. In spite of their highly simplified nature, limited mobility nonequilibrium 
surface growth models have attracted a great deal of attention [|T| primarily because of the 
following three reasons: (1) these are the first growth models which were demonstrated to 
lie outside the generic universality class; (2) in spite of their highly simplified nature, the 
various growth exponents and growth morphologies in these models seem to agree well with 
those in full diffusion Arrhenius hopping MBE growth models (where all atoms at the growth 
front are allowed to hop continuously according to temperature dependent hopping rates) 
in the low temperature kinetically rough regime; (3) obtaining a coarse-grained continuum 
description of the discrete low mobility growth models (i.e. writing down a continuum dy- 
namical growth equation corresponding to the cellular automata discrete rules defining the 
low mobility growth models) has turned out to be extremely difficult in spite of the decep- 
tive simplicity of the models. The innovation and the new technique we introduce in the 
study of these extremely well-studied limited mobility growth models is the use of the noise 
reduction technique in the direct numerical simulation of the atomistic growth rules. Our 
conclusion, based on extensive numerical simulations in both 1+1 and 2+1 dimensions (plus 
some limited simulations in higher dimensions) is that the asymptotic universality class of 
the various limited mobility growth models is a surprisingly subtle issue with many of the 
earlier findings (including some from our group) are incorrect due to pathologically slow 
crossover and extremely long transient effects which typically distract from ascertaining the 
true asymptotic universality class of these models. 

One strikingly novel feature of our results is the apparent dependence of the asymptotic 
universality class on the system dimensionality (li), not merely in the sense of well-known 
hyperscaling relations, but in the fact that the applicable hyperscaling relation itself (con- 
necting the dynamical exponent z with the roughness exponent a for example) for a specific 
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limited mobility growth model may depend on d. We believe that such dimension-dependent 
nonuniversal universality, where the hyperscaling relation for a given model changes with d 
(and which, to the best of our knowledge, has no known analog in equilibrium critical phe- 
nomena where a given hamiltonian or free energy functional, e.g. the 0^-model, is character- 
ized by a unique hyperscaling relation which does not change with the specific value of d), is 
a characteristic feature of nonequilibrium processes, and may transcend the specific models 
and phenomena being studied here. One implication of this peculiar dimension-dependent 
universality is that the universality class concept, which is one of the most important con- 
ceptual foundations of modern critical phenomena, may be of rather limited validity and 
usefulness in nonequilibrium phenomena since the model by itself does not specify the uni- 
versality class — it depends on the model (i.e. the specific set of discrete dynamic growth 
rules in our growth model) and also on the dimensionality. 

We study the limited mobility nonequilibrium growth models introduced by Wolf and 
Villain (WV) 0] and by Das Sarma and Tamborenea (DT) [Q, both of which have already 
been studied extensively in the literature and discussed rather elaborately in recent reviews 
^ . We also study two simple variants of these models, which we call asymmetric DT (ADT) 
and symmetric WV (SWV) models. The growth rules for these models are described in Fig. 
|I] for d = 1 + 1 — the rules for higher dimensions involve straightforward generalization 
of the rules depicted in Fig. |I[ In all four models atoms are deposited randomly on a 
[d — l)-dimensional square lattice substrate (which is flat initially) obeying the SOS con- 
straint — a deposited atom is then allowed to diffuse or hop instantaneously at incidence (all 
limited mobility models are also by definition, instantaneous relaxation models) to its final 
incorporation site (after which the incorporated atom never moves again) according to the 
bonding configurations of the deposition site. In DT(WV) model the deposited atom tries to 
increase (maximize) its nearest-neighbor bonding configuration by moving to the incorpora- 
tion site. If the deposition and all possible incorporation sites have the same nearest-neighbor 
bonding configurations, the incident atom does not move and stays at the site of deposition. 
The diffusion length (/), the distance over which the deposited atom is allowed to search to 
find its incorporation site, is a parameter of the model, and the simulation results shown in 
this paper use I = 1 (with the length unit being the lattice constant throughout this paper), 
i.e. nearest-neighbor diffusion only. We have, however, verified that I is an irrelevant vari- 
able, and all our conclusions in this paper are independent of the precise value of / as long as 
/ -C L, where L is the linear size of the substrate, although finite size and crossover effects 
may be strongly (and nontrivially) dependent on the value of l/L. The DT model has the 
additional diffusion constraint that only deposited atoms with no lateral nearest-neighbor 
bonds (i.e. the random incidence site has no occupied nearest-neighbor in the same "layer") 
are allowed to move whereas in the WV model all deposited atoms may move provided the 
local coordination number is maximized. The ADT and the SWV models are intermediate 
in their dynamics with respect to DT and WV models in the sense that ADT maximizes 
the local coordination while still allowing only deposited atoms with no nearest-neighbor 
lateral bonds to diffuse whereas the SWV model only increases the local coordination while 
allowing all deposited atoms to move (provided they can increase their local coordination 
or nearest neighbor bonding configuration). In all the models, the incident atom is allowed 
to move randomly to the incorporation site if several possible incorporation sites satisfy the 
growth rules. 
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All our simulations utilize the noise reduction technique introduced by two of us in an 
earlier paper The noise reduction technique allows only a fraction of the successful 
hits (i.e. the incident atoms which satisfy the specific growth rules of the model) to be 
executed — the noise reduction parameter m defines the number of successful hits needed 
before incorporation is allowed (m = 1 is the original DT or WV model). It is well-known 
that the noise reduction technique, which has been extremely successful in clarifying the 
universality class of various non-SOS growth models (e.g. Eden model, DLA) with severe 
crossover problems, is very effective in obtaining the universality class of growth models. In 
particular, the noise reduction technique (m > 1) effectively suppresses the severe correction 
to scaling induced by the strong stochasticity in limited mobility growth models, and thus 
makes it easier to determine the asymptotic universality class. We emphasize that the noise 
reduction technique is absolutely crucial in obtaining the results and conclusions presented 
in this paper, and its introduction in the simulations is the key feature which enables us to 
determine the universality class(es) of the limited mobility growth models we study in this 
paper. 

Before presenting our simulation results we briefly describe the continuum dynamical 
growth equation approach we use in discussing the universality class of various models. 
Denoting the dynamical height fluctuation variable as h, and the lateral coordinate (along 
the substrate) as x and the growth "time" as t (note that the growth time is defined entirely 
by the average deposition rate since we neglect bulk vacancies, overhangs, and evaporation 
in our SOS model) the most general leading order growth equation for our problem can be 
written as 

dh 

= u^V^h - u.V'h + + r/, (1) 

ot 

where h = /i(x, t), and r] is the spatio-temporal Gaussian white noise associated with the 
incident beam shot noise. The symmetry-allowed fourth order nonlinear term V(V/i)^ has 
been left out of Eq. (|^) since, upon renormalization, it generates the linear V^/i term already 
included in Eq. (|1]). (The constant flux term associated with the average growth has been 
left out of Eq. (|l|) since the height fluctuation variable h is defined with respect to the 
average interface position.) When 1/2 7^ 0, the asymptotic growth universality class is the so- 
called EW universality (z/2 > 0) or the "unstable" or mounded growth universality (1^2 < 0) 
— the fourth order terms in Eq. (|I]) are then all irrelevant although they may be important 
in controlling the non- asymptotic transient regime which may last for a long time. When 
z/2 = (as well as the corresponding fourth order nonlinear term V(V/i)^ being absent), 
the asymptotic growth universality class is the MBE growth universality determined by 
the V^(V/i)^ term with the initial transient regime (which could be extremely long-lived 
depending on the ratio Z/4/A22) controlled by the irrelevant fourth order linear term. When 
both z/2 = and A22 = 0, the growth universality class is the MH universality defined by 
the z/4V'^/i term. The various growth exponents for these universality classes are well-known 
and can be found in the literature. 

Determining the asymptotic growth universality class of a discrete nonequilibrium growth 
model may be difficult due to extremely long-lived transients which would mask the crossover 
behavior. This is obviously a more acute problem when any two (or all three) of the growth 
terms (z/2, A22, ^4 7^ 0) are nonzero in Eq. (|l|) although complications could also arise from 
higher order terms (i.e. higher than 4th order) neglected in Eq. ([l|). It is reasonable to 



4 



assume that in the generic situation, i.e. in the absence of any compelhng symmetry or 
conservation law induced constraints, the fourth order nonhnear equation, Eq. (|l]), should 
suffice to determine the asymptotic universality class of a given growth model, and there 
is no need to consider a higher-order (e.g. the sixth order) dynamical equation. This 
is particularly true since simple power counting considerations indicate that many of the 
higher order nonlinearities produce the V2 oi the A22 term upon renormalization, and most 
of the higher order terms are simply irrelevant. Thus it is quite possible that even in the 
extremely unlikely situation that all the terms in Eq. (|I|) vanish, i.e. V21 ^221 = 0, for some 
pathological reasons Eq. (|ip may still remain a valid generic description of the asymptotic 
growth universality class, because higher order nonlinearities neglected in Eq. (|1]) give rise 
to the terms V^/i and/or V^(V/i)^ in the growth dynamics. 

We have determined the asymptotic universality class from our discrete stochastic sim- 
ulations by calculating the growth exponents characterizing the height-height correlation 
functions as well as by measuring the surface current on a tilted substrate as proposed in 
ref . . We emphasize that the noise reduction technique is crucial in enabling us to conclu- 
sively determine the asymptotic universality class of the four (DT, WV, ADT, SWV) discrete 
growth models we study in contrast to earlier simulational studies which have been severely 
hampered by slow crossover, long transient, and strong correction to scaling problems. Our 
conclusions about the growth universality classes of the four models, as summarized in Table 
1, are based on consistent results from simultaneous measurements of growth exponents and 
surface currents on tilted substrates. 

The result md = 1 + 1 for the DT and the WV model are already known in the literature 
0J^J^,^] from extensive earlier one dimensional simulations of these two well-studied models. 
The one dimensional DT (WV) model belongs to MBE (EW) asymptotic universality as we 
also confirm decisively in our current simulations. The fact that the d = 1 + 1 DT model has 
z/2 = exactly follows from a hidden dynamical symmetry [0 in the DT growth rules which 
are symmetric, i.e., the DT diffusion rules do not have any preference between two- and 
three- bonded incorporation sites. This also makes the current on a tilted substrate exactly 
vanish in the d = 1 + 1 DT model, confirming that it is generically in the MBE universality 
class [z/2 = in Eq. (|1|)]. Similarly, earlier tilted substrate current measurements as well 
as direct simulational exponent measurements for the d = 1 + 1 WV model establish it to 
be in the EW universality class (i.e. z/2 > 0) although the WV model shows very similar 
scaling behavior to the DT model for a very long transient regime since \v2\ ~ and z/4, 
A22 are rather large in the WV model. Our noise reduced simulations of the WV model 
verify rather strikingly that this model belongs to EW universality in 1 + 1 dimensions. The 
suppression of crossover and correction to scaling effect in the 1 + 1- dimensional DT and 
WV models was our original motivation for introducing the noise reduction technique in 
this context. 

The kinetically rough surface growth is traditionally analyzed in terms of the dynamic 
scaling hypothesis where the surface width (the root mean square fluctuation in the surface 
height) or more generally, the height-height correlation function shows generically scale 
invariant power law scaling behavior, with critical exponents a, /?, z = a/f3, given by 
W{L,t) ~ L°'f{^{t)/L) where W is the dynamical surface width at growth time t for a 
substrate of lateral size L and ^ (t) is the lateral dynamical correlation length for the speciflc 
growth process with ^{t) ~ t^^^ . The scaling function f{x) behaves as f{x ^ 1) ~ 1 and 
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/(x < 1) ~ a;" so that W{L, t ^ oo) = W,{L) ~ L° and W{L, t < L^^) ~ tl^ with f3 = a/z, 
where Ws(ij) is the saturated steady-state long time surface width and W{L <^ Hf)) is the 
pre-steady-state dynamical surface width. The set of critical exponents a, /3, and z = a/P 
define the growth universality class, and in principle, can be determined for a particular 
continuum growth equation. For the sake of completeness (see the WV growth results 
presented below) we point out that for mounding instability with slope selection, when the 
surface morphology evolves into a regular mounded pattern with the sides of the mounds 
having constant slopes, one gets a = 1, and only one exponent {(3 oi z) (3 = 1/z then 
defines the growth pattern. In Table 1 and in Figs. ^ and ^ we present our calculated 
critical exponents for the DT and WV model along with the associated continuum equation 
descriptions. 

Our most surprising findings are presented in Fig. ^ (for the DT model) and Fig. ^ (WV 
model) where our results for the 2 + 1- dimensional DT and WV models are depicted. We 
find that the DT (WV) model in 2 + 1 dimensions belongs to the EW (unstable) dynamic 
universality in contrast to the 1 + 1 dimensional universality class of these models. This is the 
important new result presented in this paper which disagrees with the earlier conclusions 
in the literature. The determination of the asymptotic universality class of the 2 + 1- 
dimensional DT and WV model is the main new result being presented in this paper. 

We first discuss the 2 + 1 dimensional DT results which are in some sense less surprising 
than the corresponding WV results because EW universality is the generic SOS universality 
class. Our measurement of current on tilted substrates in the 2 + 1 dimensional noise reduced 
DT simulations always exhibit a small but finite downhill current indicating the presence of 
a small z/2 > term in Eq. (|lD. This should be contrasted with the corresponding 1 + 1- 
dimensional DT results where a simple symmetry argument as well as extensive numerical 
simulations definitively establish the absence of the V^/i term in the d = 1 + 1 DT growth 
equation, i.e. 2^2 = 0(7^ 0) in the d = 1 + 1(2 + 1) DT model. The asymptotic universality 
class of 1 + 1- dimensional DT model is now well-established to be given by the following 
continuum growth equation, which does not have the generic 1^2^ term of Eq. (0): 

where the "higher order" terms of the form ^(f^)^™ with 2m = n = 4,6,8,..., etc. are 
marginally relevant in d = 1 + 1 (i.e. simple power counting reveals them as having the same 
anomalous dimensions as the nonlinear fourth order A22 term). Our tilted substrate current 
measurement in c? = 2 + 1 DT model shows the existence of a small slope dependent downhill 
current of appoximate magnitude ~ 10~^ whereas the corresponding d = 1 + 1 DT model 
has a current ~ 10~^ (of random sign) which is indistinguishable from the background noise 
effect. We therefore conclude that the 2+1 dimensional DT model, describing nonequilibrium 
growth on physical surfaces and interfaces, belongs to the generic EW universality class, and 
has the following coarse-grained continuum description: 

dh 

= y,V^h-v,V'h+ J2 A2(2n)V2(V/i)2" + r/. (3) 

n=l,2,3,... 

Our finding that z/2 is very small, but nonzero, in Eq. (^ for 2 + 1- dimensional DT growth 
should not come as a big surprise (except, of course, for the fact that it has not earlier been 
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discovered in the literature including our own earlier work on the DT model) because 
the vanishing of z/2 in the 1 + 1- dimensional DT model arises from a rather peculiar kinetic- 
topological symmetry of the DT model, which applies only in one dimension and cannot be 
generalized to two dimensional surfaces. In the absence of a compelling symmetry argument 
manifestly making z/2 = in the growth equation, one should expect its generic presence in 
the 2 + 1- dimensional DT model although the extreme quantitative smallness of z/2 has made 
it difficult so far to establish its finiteness in simulations. Our finding that there is a small 
downhill current on tilted substrates in the 2 + 1- dimensional DT model and that the critical 
exponents of 2 + 1- dimensional DT growth (Fig. 0) are consistent with EW universality 
class (and not particularly consistent with the MBE universality defined by z/2 = in Eq. 
(HI)) leads to the conclusion that 2 + 1- dimensional DT growth is in the EW universality 
class (z/2 7^ in Eq. {^) and 1 + 1- dimensional DT growth is in the MBE universality class 
(z/2 = 0). We have also carried out DT simulations in (unphysical) 3 + 1 dimensions finding 
very good agreement with EW universality properties. 

Our results for the 2 + 1- dimensional WV model (Fig. |^) are very dramatic and com- 
pletely unanticipated. We find that the 2 + 1- dimensional WV model leads to spectacular 
quasiregular mounded morphology indicating unstable epitaxial growth (whereas the corre- 
sponding 1 + 1- dimensional WV growth is asymptotically stable and fiat, belonging to the 
EW universality class). Thus the usual critical exponents {a, (3, z = a/P) are not particu- 
larly meaningful for 2 + 1- dimensional WV growth (although they can still be defined in the 
simulation results — the exponents, however, provide a misleading picture since the growth 
front, instead of being statistically scale invariant as it should be for kinetic surface rough- 
ening exhibiting power laws controlled by critical exponents, has a quasiregular mounded 
pattern). 

The tilted substrate current measurement in the 2 + 1- dimensional WV growth yields 
another curious surprise. It turns out that the local current is uphill (along 111 plane) 
or donwhill (along 100) depending on how (i.e. along which direction) one decides to tilt. 
The fact that the tilted substrate current could depend sensitively (i.e. stabilizing in some 
directions and destabilizing in other directions) on the tilt direction has not earlier been 
reported in the literature where most reported current measurements are carried out in 
1 + 1 dimensions (where, of course, this problem cannot arise) with the hope or expectation 
(proven to be false in this paper) that an accurate determination of the universality class of 
a growth model in 1 + 1 dimensions will automatically give us the same universality class in 
2+1 dimensions — WV model is in the EW universality class in 1+1 dimensions and unstable 
(mounded morphology) in 2 + 1 dimensions! Thus the tilted substrate current measurement, 
while being capable of providing the correct universality class in 1 + 1 dimensions, may very 
well lead to misleading and wrong conclusions in higher dimensions where the current on 
tilted substrates is explicitly direction dependent and is not uniquely defined. We will publish 
quantitative details on the WV mounding phenomenon elsewhere — here we point out that 
(1) the underlying mechanism for the WV mounded morphology is related to surface cluster- 
edge (or kink) diffusion induced mounding recently discussed in the literature |0, and (2) 
this WV mounding phenomenon, arising as it does from kinetic-topological aspect of surface 
diffusion, leads to very strong instabilities in (unphysical) dimensions higher than 2+1, where 
early work reported |jlO[ unexplained strong mound formation in 3 + 1 and 4 + 1 dimensional 
WV growth. We have carried out WV simulations in 3 + 1 dimensions, finding very strong 
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mounding even without any noise reduction consistent with earher findings |T0 



Finally we consider the two intermediate models, ADT and SWV, for the sake of com- 
pleteness. The d=2+l growth morphologies in these two models are shown in Figs. ^ (ADT) 
and 1^ (SWV). Without presenting the actual numerical results for the critical exponents for 
these two models, we just mention that the ADT and SWV models in 2+1 dimensions be- 
have qualitatively similar to the WV model in 2 + 1 dimensions with fairly strong mounding 
under noise reduction although the morphological details in the two models differ somewhat 
with the SWV morphology being similar to the pyramidal structures of the WV morphology 
and the ADT morphology having flat top mounds with very deep and narrow grooves. Thus, 
ADT/SWV/WV all have unstable growth in 2 + 1 (or higher) dimensions with quasiregular 
mounded morphology whereas DT in 2 + 1 (or higher) dimensions is in the EW universality 
with stable and smooth growth morphology. All four models have kinetically rough statisti- 
cally scale invariant growth in 1 + 1 dimensions with the WV (DT) model belonging to EW 
(MBE) universality class. 

We conclude by stating that we have found the universality class concept to be of limited 
usefulness in conserved discrete limited mobility nonequilibrium surface growth models. The 
same growth rules deflning a particular model (e.g. WV or DT) may belong to different 
universality classes in different dimensionalities (not in the sense of super-universality, but 
in a more fundamental nontrivial sense as if an equilibrium model, which is in the Ising 
universality class in two dimensions, behaves as an x-y model in three dimensions — a 
patently absurd notion). In addition, rather minor changes in local growth rules could lead 
to dramatic differences in the resulting growth morphology or the universality class — DT 
and WV have very similar local growth rules, but their morphologies, smooth (DT) and 
mounded (WV), and universality class (EW for DT and unstable for WV) are strikingly 
quahtatively different. We also flnd that measuring surface current on tilted substrates 0, 
while being a potentially useful technique for discussing the 1 + 1 dimensional universality 
class, may not work in 2+1 dimensions and may produce misleading or conflicting conclusions 
depending on the precise direction of the surface current. We have assumed throughout that 
the noise reduction technique, which is absolutely crucial in our obtaining the asymptotic 
universality classes of various models we study, does not modify the universality class of a 
growth model (and only suppresses transient and correction-to-scaling effects by reducing 
the effective noise strength). This belief is based on extensive earlier analysis of the noise 
reduction technique in the literature, which, in general, is thought not to affect the growth 
universality class. 

This work is supported by the NSF-MRSEC and US-ONR. 
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FIGURES 



FIG. 1. Schematic configuartions defining the growth rules for the DT and WV models in 
d = 1 + 1. For ADT model, adatom "g" in the top (DT) plot will only hop to the left where it 
will have 3 nearest neighbors. For SWV model, adatom "g" in the bottom (WV) plot will choose 
between the left and right neighbors with equal probability. 

FIG. 2. (a) The interface width W versus time t in 2+1 DT model with substrate size 
L = 1000 X 1000 and noise reduction factor m = 1, 3, 10, and 15 from top to bottom. Solid lines 
are best power law fit which yield the growth exponent /? that decreases as m increases, (b) A 
typical morphology of a noise reduced DT model (m = 5) from a 1000 x 1000 substrate (only a 
section of 200 x 200 is shown here) at 400 ML. 

FIG. 3. (a) The interface width W versus time t in 2+1 WV model with substrate size 
L = 100 X 100 and noise reduction factor m = 5. (b) A typical morphology of a noise reduced WV 
model (m = 5) from a 500 x 500 substrate (only a section of 200 x 200 is shown here) at 10^ ML. 

FIG. 4. A typical morphology of a noise reduced ADT model (m = 5) from a 100 x 100 
substrate at 10^ ML. 

FIG. 5. A typical morphology of a noise reduced SWV model (m = 5) from a 100 x 100 
substrate at 10^ ML. 
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TABLES 



model 


properties 


d=l+l 


d=2+l 


DT 


exponents 


a = l, (3 = 1/3, z = 3 


a = O(log), P = O(log), z = 2 


morphology 


kinetically rough 


flat and very smooth 


particle current 


J = 0(~ ±10-*^) 


J ~ -10-4 


WV 


exponents 


a = 1/2, /3 = 1/4, z = 2 


a = l, p= 1/3, z = 4 


morphology 


kinetically rough 


mounds with selected slope 


particle current 


J ~ -10-3 


Jioo - -10-2; Jni ~ +10-2 



TABLE I. Summarized of our results for the DT and WV models 
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a 



DT model 



a 



WV model 



